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Abstract 

A fast method to determine the crystallographic axes of segmented true-coaxial 
high-purity germanium detectors is presented. It is based on the analysis of segment- 
occupancy patterns obtained by irradiation with radioactive sources. The measured 
patterns are compared to predictions for different axes orientations. The predictions 
require a simulation of the trajectories of the charge carriers taking the transverse 
anisotropy of their drift into account. 
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1 Introduction 



High purity germanium detectors, HPGeDs, are used in a wide variety of 
applications in particle and nuclear physics [Tp|3|H] . Segmented detectors can 
be used to disentangle event topologies [5|l6|7] . The detectors considered here 
are segmented cylindrical true-coaxial detectors. The segmentation in azimuth 
angle [IJ is the decisive feature for the proposed method. 
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The orientation of the crystallographic axes of a germanium detector is impor- 
tant for many analyses where pulse shapes are used. The difference between 
the charge-carrier mobilities along the crystallographic axes (110) and (100) 
is significant. It is called longitudinal anisotropy. Experimental data on hole 
mobility exist for the (111) and (100) axes [8]. Relevant here is the plane con- 
taining the (110) and (100) axes. Appropriate calculations [9] imply that the 
difference of the velocities of the holes along the (110) and (100) axes is about 
10 % for realistic values of the electric field. The different charge carrier veloc- 
ities result in different rise times of pulses originating at equal distances from 
the electrode. The differences in rise times between pulses along the (110) and 
(100) is typically around 10% PITO] depending on the relative location of the 
charge deposits to the crystal axes. In addition, the paths of charge carriers 
not drifting along the crystallographic axes are bent; this is called transverse 
anisotropy. The anisotropics have to be taken into account in so called pulse 
shape analyses [TTP^ and whenever simulated pulses [S] are compared to data. 

The crystallographic (OOl)-axis of a cylindrical germanium detector is usually 
aligned with the z axis. The position of the (110) axis, 0{iio>, is, however, a 
'priori not known after the detector is processed. Therefore the orientation of 
the (110) axis with respect to the segment boundaries has to be determined 
during the characterization of the detector. 

This paper presents a comparison of two methods to determine the axes. The 
first method is based on scanning the detector. It is widely used, well proven 
and based on data only. However, it requires that the detector is mounted in a 
special test setup with a movable source and it is time consuming. The second 
method involves the comparison to Monte Carlo simulation. Its advantages are 
that it is fast and can be performed in any configuration, even if the detector 
is part of a complex detector system. The determination of the axes is required 
to have an accuracy of about 5°. A simple application using the orientation of 
the axes is the calculation of the probability that a neutrinoless double beta 
decay event actually produces a signal in two segments of a detector, even 
though it is a so called single-site event. This requires that the knowledge 
about the axes orientation is sufficient to calculate whether a drifting charge 
cloud gets separated across segment boundaries. The charge cloud representing 
a neutrinoless double beta decay event has a diameter of about 2 mm [13], and 
thus covers about 5° at a radius of 25 mm in a cylindrical detector. 



2 Detector and test environment 

The detector used for this study was an n-type true-coaxial 18-fold segmented 
detector with a diameter of 75 mm, a height of 70 mm and an inner bore of 
10 mm. It was produced by Canberra, France. The segmentation scheme is 
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three in z and six in 0. The layout is depicted in Fig. [T] The first such detector 
was characterized in detail [lO] including the determination of the orientation 
of the crystallographic axes. The second detector of the series, used for this 
study, was previously used to test the performance of such a device while 
submerged in a cryogenic liquid [H] and the temperature dependence of pulse 
lengths [T5] . 



^180- 




a) b) 

Fig. 1. (a) Coordinate system with the x axis chosen to be on the 4-15 segment 
boundary. Segment boundaries are indicated as dashed lines, coordinate 
dotted lines. Numbers in parenthesis are segment numbers as in (b). (b) Detector 
segment layout in (0,2:) coordinates. 

The density of electrically active impurities, pim-p-, was given by the manufac- 
turer as 0.35 ■ 10^°cm~^ at the top and 0.55 ■ 10^°cm~^ at the bottom of the 
detector. The change in impurities is assumed to be linear with height, z. The 
operational voltage was 2000 V. At the time of the measurements, the detector 
drew a leakage current of up to 100 nA, all through segment 9. It was, however, 
fully functional. The only effect was that the energy resolutions in segment 9 
and the core had deteriorated. 

For all measurements, the detector was mounted inside the vacuum cryostat 
Kl [15]. The detector was scanned in at 2 = with a ^^^Eu source placed 
at r=9.6cm. It was also irradiated from the top (z = 18.5 cm, r=Ocm) with 
a ^^^Th and a ^^Co source. The ^^'^Th source was also placed at the side 
(2; = —2.3 cm, r = 17.6 cm, = 155°), in front of segment 16. 



3 Detector scan 

A well established way to find the crystallographic axes is to perform an 
azimuthal scan, where a source with a low energy gamma line, like ^^^Eu with 
its 122 keV line, is moved around the detector in small steps in 0. Low energy 
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gamma rays predominantly deposit their energy close to the surface. The rise 
times of the resulting pulses vary with 0, reaching a maximum for the drift 
along the (110) axes. 

A scan was performed with a collimated 40 kBq ^^^Eu source with a la beam 
spot diameter of about 5 mm and a step size of 10° in 0. The events were 
selected using an energy window of 10 keV. The rise time of a given pulse was 
determined by fitting a simulated pulse to the measured one. This is described 
in detail elsewhere piTS] . This method uses all the available information about 
a measured pulse while suppressing the influence of noise. For every position, 
the average rise time, tr, was computed for all pulses, for which x^/ndof < 1.5 
for the fit with the simulated pulse. Figure [2] shows the dependence of tj. on 
the angle of the source position, . The data were fitted with the function 



where A, B and 0(iio) are the free parameters. The coordinate system was 
chosen such that the 4-15 segment boundary was at 0°. 

A study of the systematic uncertainties was performed. The energy window, 
in which events were selected, was changed to 5keV and 20keV. In addition, 
the used for the pulse selection was modified from 1.5 to 1.2 and 1.8. 
In both systematic uncertainty of 0.5° was deduced. The dominating 

systematic uncertainty is connected to the placement of the detector within the 
cryostat. This could only be controlled to 3°. The total systematic uncertainty 
was determined by adding the individual contributions in quadrature. The 
statistical uncertainty, as determined by the fit, is 0.4°. The final result is: 



4 Occupancy method 

The drift of the charge carriers within a germanium detector only follows 
radial lines along the crystallographic axes. The transverse anisotropy causes 
bent trajectories elsewhere in the crystal [9l Fig. 3]. As a result, the effective 
segment boundaries do not coincide with the geometrical boundaries and the 
segment electrodes collect charge carriers from volumina of different sizes. This 
reflects in the occupancies of the segments when the detector is irradiated. 




(1) 



0(110) 



0.2° ±0.4°(stat.) ±3.1°(syst) 



(2) 
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Fig. 2. Dependence of rise time on the azimuth angle (p of the source position. The 
points are the data with statistical uncertainties. If the uncertainty bars are not 
visible, they are smaller than the points. The line represents the fit described in the 
text. The vertical lines indicate the segment boundaries; the numbers are segment 
numbers. 

4-1 Measured Occupancies 



A 28kBq ^^^Th and a 40kBq ^''Co source were used. The occupancies were 
measured for the 0.58 MeV and 2.61 MeV lines from ^o^Tl and the 1.17MeV 
and 1.33 MeV lines from ^''Co. Only events, in which one segment registered 
the same energy, Es, as the core within 50keV, and the energy seen in any 
other segment was less than 150 keV were used. As the separately measured 
background spectra did not show lines associated to ^°Co, the background 
was not subtracted for the analysis of these lines. The background did show 
lines associated to ^^^Tl decays. Therefore the background was subtracted for 
the analysis of the corresponding spectra. Figure [3]^a) shows the spectrum 
measured in segment 3 for the irradiation from top and the corresponding 
normalized background spectrum. The background contribution for this rela- 
tively low-energy line is maximal at the bottom of the detector. The resulting 
segment-spectra were fitted using the log-likelihood method. The fit function 
was a combination of a Gauss function of width a and a sigmoid: 

f{E,) = ^^e ^+5 + ^ , (3) 



where E^ is the energy of the line under study and A, B, C are the free pa- 
rameters. The sigmoid represents the background shape. It drops for energies 
Eg > E^. The occupancy, Di, of each segment i was taken as the fitted number 
of events under the peak, = A - binwidth"^ 

Figure [3](b) shows the result of the fit to the background-subtracted spectrum 
around 0.58 MeV for the irradiation with ^ostI from the top. The contributions 
from the terms associated to the parameters A, B and C are indicated. The 
area in the peak associated to the parameter A is well defined and what enters 
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Energy [MeV] Energy [MeV] 



Fig. 3. (a) Measured signal and background spectrum in segment 3 around 0.58 MeV 
and (b) the corresponding fit to the background-subtracted measured spectrum with 
the function given in Equation [3l 
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Fig. 4. Energy spectra at 1.33 MeV as seen by the 18 segments for the irradiation 
with ^'^Co from the top. Also shown are the results of the fits with the function 
given in Equation [3l 

the analysis. Figure H] shows the spectra of all 18 segments for the irradiation 
with ^"^Co from the top Also shown are the fit results. 

Figure [5]^a) shows the resulting occupancies for the middle layer. A clear struc- 
ture is visible. It is almost mirrored with respect to the boundary between 
segments 13 and 6. This agrees with the naive expectation, because the (110) 
axis is almost aligned with this boundary. 



Segment 9 in the top layer has a significantly larger resolution and its occupancy 
is only used in the analysis of data obtained by irradiation from the top. 
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Fig. 5. Measured occupancies (a) for the 1.33 MeV line extracted for the middle 
layer for the irradiation from the top and (b) for the 2.61 MeV line extracted for 
the middle layer for the irradiation from the side. The numbers denote the segment 
numbers. 

Figure [5](b) shows the measured occupancies for the 2.61 MeV line in the mid- 
dle layer for the irradiation from the side. The situation is more complicated in 
this case and there is no naive expectation for the influence of the transverse 
anisotropy. The occupancies were extracted for all layers and all energies and 
for both irradiation from the top and the side. 



4.2 Expected Occupancies 



Events for the experimental setup as described in Section |2] were simulated 
with the MaGe [16] package and trajectories were computed with the pulse 
shape simulation package [9]. The package provides calculations of the electri- 
cal and weighting fields and uses experimental input and model based calcula- 
tions to provide the mobilities, and thus the velocities, of holes and electrons 
at any given point in the detector. The orientation of the (110) axis, 0(no>' 
was an input parameter to the simulation of the drift. The endpoints of the 
resulting trajectories were used to assign each energy deposit of an event to 
a segment. The energies were summed for each segment. The event selection 
was performed as for the measurements. 

The resulting expected occupancies for (/>^'™o^ = —20° and (t>^(^{^o) ~ 
shown in Fig. [6] for the irradiation from the top and the side. The occupancies 
depicted in Fig. [6] were computed for the middle layer, for which an average 
Pimp = 0.45 ■ 10^°cm~^ was assumed. For the irradiation from the top, the 
result for the 1.33 MeV line of ^°Co is shown. For the irradiation from the 
side, the 2.61 MeV line of ^os^l was chosen. 

For the irradiation from the top, the simulation of the middle layer predicts 
occupancy patterns which are almost identical for segments 15-14-13 and 
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Fig. 6. Expected occupancy patterns in the middle layer for the irradiation from 



the top/side for (a)/(c) = -20° and for (b)/(d) ^f^^^^ 



30°. 



6-5-4. This reflects the 180° degeneracy of the axis orientation. The mirror 
symmetry observed in data with respect to the 13-6 boundary disappears, if 
the (110) axis is not aligned with this boundary. 

For the irradiation from the side, a clear difference between the predictions for 

30° is visible. There are no apparent symmetries. 



'(110) 



-20° and 



'(110) 



but the relative occupancies significantly change for varying 0(iio)- 

The amplitudes of the observed patterns depend on the amount of transverse 
anisotropy in the crystal. This is influenced by the hole mobility and by pimp- 
The dependence on pimp was investigated by comparing expectations for the 
nominal detector, i.e. the one supposed to be close to the real device with a 
Pimp of 0.52 ■ 10^°cm"^, 0.45 ■ 10^°cm"^ and 0.38 • 10^°cm"^, to predictions for 
a lower impurity test detector with a pimp of 0.29 • 10^°cm~^, 0.18 • 10^°cm~^ 
and 0.06 ■ 10^°cm~^ for the top, middle and and bottom layers. Figures [^a) 
and (b) show the expectations for 



'(110) 



0°. 



The amplitude of a pattern can be defined as the depth of the first step of 
the pattern, Ai_j = Si — Sj, with Si and Sj being the predicted occupancies 
for the segments i and j located in the same layer. For the middle layer, 
^15-14 = 0.0326 is predicted for the nominal detector and Ai5_i4 = 0.0336 for 
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Fig. 7. Comparison of the normalized predicted occupancies for the 1.33 MeV Une 
and ^|™o^ = 0° for the middle, top and bottom layers of (a) the nominal crystal 
and (b) a test crystal with lower pimp- 

the test detector. This is a very small change of 0.6 % in a 20 % effect with 
respect to an occupancy of S'j = 1/6 = 0.17 for all i, which would be seen if 
there was no transverse anisotropy. 

The Figs. [Tl^a) and (b) also show predictions for the top and bottom layers for 
the two simulated devices. Shown are deviations, A^.^ = Sk — Si with k and 1 
being the relevant segment numbers in different layers, but at the same 0. The 
values of A^.; reach 0.003, which is a 2 % effect, i.e. 10% of the amplitude. 
The 15-14-13 6-5-4 degeneracy is broken for the end layers due to a drift 
component in z, which arises from pimp changing with z. The Afc_; patterns 
are basically identical for the nominal and the test device. 
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Fig. 8. Expected occupancies for (f'f^^Q) = 0° ^oi the middle layer for (a) the 
1.33 MeV line and for (b) the 2.61 MeV line for irradiation from the top. Also shown 
are the changes predicted for a 5 mm misalignment of the source along 45° or 90° . 

The sensitivity of the predictions on the exact location of the source is demon- 
strated in Fig. [8] for the case of (a) the 1.33 MeV line of ^"Co and (b) the 
2.61 MeV line of ^°^T1. Shown are the normalized predicted occupancies for a 
perfectly aligned source and the effects on the pattern for two source misalign- 
ments of 5 mm, one along = 45° and one along = 90°. The misalignments 
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cause changes in the pattern for the middle layer breaking the 15-14-13 6-5-4 
degeneracy. The effect on a single occupancy is as much as 0.005, which is a 
3 % effect on the occupancy, i.e. as much as 15 % on the amplitude. 



4-.3 Extraction of the axes orientation 



The occupancies were computed for 0™o> varying in 1° steps. The resulting 
occupancies were analyzed by computing the test statistic e, defined as 



E 



(A - S,) 



(4) 



where Di and Si denote the measured and simulated occupancies in segment i. 
The simulation was normalized to the total number of events in a given layer 
and the sum runs over all segments in this layer. 



The resulting function e(0^no>) fitted in a 20° window with a second order 
polynomial. The 0^i™o> corresponding to the minimum of the polynomial, emin, 
was taken as the result of the procedure, ^"'^'^.s^ Figure [9](a) shows e(0^'™Q^) 
plus the result of the polynomial fit for the case of the 1.33 MeV line in the 
middle layer for irradiation from the top. The shape of the pattern observed 
in the data is best reproduced by the prediction for 0^™o) ~ 4°. Figure Mjo) 
depicts the corresponding measured and predicted occupancies. 
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Fig. 9. (a) The dependence of e on <A^™o^ for the 1.33 MeV line in the middle layer for 
the irradiation from the top. (b) The corresponding measured and expected occu- 
pancies. The numbers denote the segment numbers. The simulation was normalized 
to the observed number of events in the middle layer. 

The comparison between the predicted and measured pattern in Fig. [9] shows 
that the shape of the pattern is very well reproduced. However, the amplitude 
of the predicted pattern is larger than the amplitude of the measured pattern. 
As was shown in Section 14. 2[ this cannot be explained by the choice of pimp 
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in the simulation or by any misalignment of the source. As the hole mobility 
is not well known [8], its simulation could be the source of the discrepancy. 



(a) 



ijj 



CO 

.03 



10 



r 


■niddl< 


3 laye 


r 






















1.11.1, 














/ 































o 

03 
Q. 

O 

o 
O 



120 140 160 



180 

^<110> 



200 
[deg.. 



800 
600 
400 
200 




(b) 

SIDE, middle layer 



— Measured 

— Simulated ' 1 



15 14 



13 



Segment 



Fig. 10. (a) The dependence of e on (/"^^q) 1A\ MeV line in the middle 

layer for the irradiation from the side, (b) The corresponding measured and ex- 
pected occupancies. The numbers denote the segment numbers. The simulation was 
normalized to the observed number of events in the middle layer. 

Figure [10] shows the situation for the 2.61 MeV line in the middle layer for 
the irradiation from the side. The value of emin is significantly higher than for 
the top case depicted in Fig. [9l but the minimum is as distinct. The pattern 
observed in the data is well reproduced by the prediction for 0(iio) ~ 164°= — 
16°. 



^.^ Systematic Uncertainties 



The systematic uncertainties fall into two categories: 

(1) parameters used in the treatment of the data 

(2) imperfections of the simulation 

The first category is dominated by uncertainties arising through the definition 
of e and the fit to the resulting e-function. The test statistic e deviates from 
a classical by putting extra emphasis on bins with lower occupancies. For 
the irradiation from the top, the differences in the results are small and the 
resulting systematic uncertainty is 1°. For the irradiation from the side, the 
differences in the results are significant. This is intended. However, a study 
using different test statistics yielded a systematic uncertainty of 4°. The effect 
of a variation of the width of the fit window for the e-function was also stud- 
ied. The window was widened to 30° and narrowed to 10° and the resulting 
uncertainty is 1°. Using other functions then a second order polynomial to fit 
the e-function yielded a systematic uncertainty of 2°. 
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The total systematic uncertainty connected with e is evaluated by adding the 
individual uncertainties in quadrature. It is 2.5° and 4.6° for irradiation from 
the top and the side, respectively. 

In the second category, the influence of the combination of hits before drifting 
the charge carriers was studied. Per default, hits within a distance of 1 mm 
were clustered. That was changed to 2 mm and the uncertainty due to this was 
found to be negligible. Any diffusion of the charge cloud is not simulated; the 
effect is expected to be small compared to the influences of charge clustering. 

The main uncertainties are connected with the assumptions on the detector 
parameters and the source location. 

As shown in Section U^l a mismatch of simulated and real pimp slightly changes 
the amphtude of the pattern, but not the shape. Therefore, the absolute val- 
ues of e will change, but the result will actually not. This was confirmed by 
simulation. 

Also discussed in Section 14.21 was the influence of any misalignment of the 
source. The source position was controlled to about 5 mm. Different layers are 
affected slightly differently due to the difference in angular coverage. The sim- 
ulations show that the systematic uncertainty due to possible misalignments 
of up to 5 mm is 5° for irradiation from the top. 

The systematic uncertainty for the irradiation from the side directly reflects 
the uncertainty in the 0-position of the source with respect to the detector. 
The detector was remounted for the side measurement and its relative rotation 
within the cryostat could, as mentioned before in Section [3l only be controlled 
to 3°. In addition, the position of the source relative to the cryostat could only 
be controlled to 5°. Added in quadrature, this results in a total systematic 
uncertainty of 6°. 

For the analysis of lines which are also present in the background, possible 
changes in the background over time are another source of systematic uncer- 
tainty. Such changes can occur, if the cryostat or any other equipment in the 
lab is moved. Unfortunately, this was the case. The resulting patterns are very 
sensitive, especially for the irradiation from the side. The uncertainty was eval- 
uated conservatively by performing the complete analysis for the ^os^jj lines 
with and without background subtraction. The result is a systematic uncer- 
tainty of 1° and 7° for the top and side cases, respectively. 

In total, the systematic uncertainties for the irradiation from the top and the 
side were evaluated to be 6° and 10°, respectively. They are dominated by 
the uncertainty in the source position and, for the irradiation from the side, 
possible changes in the background. 
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4-5 Results 



The axes orientation was determined for all energies and all layers. The re- 
sults, 0™^^^^ are given in Table [TJ Also given are weights, w = Cmirn 
measurements. The statistical uncertainty of the fit to the e-function was al- 
ways negligible compared to the systematic uncertainty. 
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Irradiation from the top: 


0.58 


-11.1° 


79 


-8.6° 


52 


13.1° 


52 


1.17 


-7.7° 


100 


1.4° 


115 


-3.7° 


167 


1.33 


-7.1° 


88 


4.2° 


99 


1.1° 


117 


2.61 


0.3° 


38 


-6.2° 


23 


-6.3° 


44 


Irradiation from the side: 


0.58 


-21.9° 


33 


-14.2° 


52 


-23.4° 


63 


2.61 


-18.6° 


19 


-15.6° 


43 


-15.9° 


24 



Table 1 

Extracted (llO)-axis positions, cj/^^^, for different gamma lines and layers. Also 
listed are the weights, w, as defined in the text. 



For the irradiation from the top, the results for the ^^^Tl lines show a spread 
which is slightly larger than expected from the evaluation of the systematic 
uncertainties. This could be due to the effect of statistical fluctuations in 
individual bins. An ensemble test was used to determine the statistical uncer- 
tainties. The test was performed by randomly varying the content of each bin 
according to its statistical uncertainty and extracting 0™'=^'^ for the modified 
occupancy pattern. The resulting 0™^^*^ distribution was fitted with a Gaus- 
sian function and the result used to determine the statistical uncertainty. The 
largest uncertainties were observed for the 2.61 MeV line for the irradiation 
from the top. They were 4° for all layers. 

The statistical and systematic uncertainties are different for different energies 
and layers. Therefore, weighted averages were computed: 



J2w 



'{110) ~ ^^77 ' (^) 



where the sums run over a combination of layers and line energies. The results 
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for the ^°Co and ^os^j lines in the three layers for irradiation from the top are: 
= -1.8° ± l°(stat) ± 6°(syst) , (6) 
(i\m = -3-6° ± 4°(stat) ± 6°(syst) . (7) 

Figure [11] shows the measured occupancies together with the best predictions 
for the middle and bottom layer for the ^°Co irradiation. The high statistics 
^°Co distributions show that the predictions reproduce the shapes of all four 
measured patterns well, but that the predicted amplitudes are consistently too 
large. The measured amplitudes are about 40 % smaller than the predicted 
ones. As discussed previously, this cannot be explained by an impurity level 
assumed to be too high in the simulation. The reproduction of the shapes of 
the measured distribution also indicates that the alignment of the source was 
well within the 5 mm uncertainty. 

The simulation uses the hole mobilities as calculated by the pulse shape simula- 
tion package [9] using measured input parameters [8] . There is a non- negligible 
uncertainty in these calculations and the input parameters. The observation 
of less pronounced patterns than predicted could eventually be used to test 
the assumptions about hole mobilities. A simulation confirmed that the dis- 
crepancy between the predicted and observed amplitudes could be explained 
by an underestimate of the ratio of the mobilities along (111) and (100). The 
trajectories were recalculated for mobilities, for which the component along 
the (111) axis was globally, i.e. independent of the electric field, increased by 
10%. This resulted in amplitudes about the size observed. However, as the 
calculation of mobilities along the (110) axis and between the axes depend on 
several approximations, this will need further systematic studies. 

Figure [12] shows the measured occupancies together with the best predictions 
for the middle and bottom layer for the "^^^Tl irradiation from the top. The 
statistics is much lower than for the irradiation with ^°Co, such that single 
bins can be affected by significant statistical fluctuations. The 2.61 MeV line 
does not have a large probability to get fully absorbed in the volume of a 
single segment. The bottom layer has fewer events by design, because of the 
absorption in the crystal. It should be noted that the spread of the single 
measurements is still reasonable, even though single bins are fluctuating. 

The individual values obtained for the irradiation from the side deviate con- 
sistently by more than 10° from the expectation. The ensemble test for the 
2.6 MeV line, as discussed previously in this chapter, yielded the statistical 
uncertainty of 4° and 6° for the middle and bottom layers, respectively. Seg- 
ment 9 was not used in the measurements using the top layer. Even though this 
segment was directly facing the source, the resulting measurements are still 
consistent with the results for the other layers. The top layer was, however. 
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Fig. 11. Measured occupancies together with the best predictions for the irradiation 
with ^'^Co from the top for (a) 1.17 MeV, middle layer (b) 1.33 MeV, middle layer 
(c) 1.17MeV, bottom layer (d) 1.33 MeV, bottom layer. 

excluded from the averaging. The overall result is: 



/Tl,side 
'(110) 



-17.9° ± 5°(stat) ± 10°(syst) . 



(8) 



This result is barely consistent with the one obtained with the scan. Figure 
depicts the measured occupancies together with the best predictions for the 
middle and bottom layer for the ^os^l irradiation from the side. An irradiation 
directly from the side is a priori not favorable for this method. However, the 
patterns are reproduced quite well, indicating that probably the misalignment 
caused by the remounting of the detector and the changes in the background 
due to the movement of the equipment were larger than hoped for. The usage of 
this data for the analysis presented here was originally not planned. Therefore 
the systematic uncertainties relevant for this analysis were not really well 
controlled. With some care, they certainly could be reduced significantly. 

The overall results confirm the ability to extract the axes orientation from 
the occupancy patterns in a 0-segmented detector, even if the detector is not 
in a perfect condition. The precision depends on the care taken to avoid sys- 
tematic uncertainties. The source position has to be well controlled under all 
circumstances. A symmetric configuration reduces the systematic uncertain- 



15 



(a) 



(b) 



r 

o 

Q. 1000: 

o 
o 

o 



500- 



TOP, middle layer 

r" 



1_ 



r ' 

— Measured 

— Simulated 



15 14 



13 



_5 4 

Segment 



o 

CO 
Q. 

o 
o 

O 



400- 



TOP, bottom layer 



300: 



200- 



100- 



Measured 
Simulated 



18 17 16 



1 



(c) 



o 

CO 
Q. 

O 

o 
O 



o 

CO 
Q. 

o 
o 

O 



Segment 



400 
300 
200 
100 


200 
150 
100 
50 



TOP, middle layer 

— 



—I I 

— Measured 

— Simulated 



15 14 



13 



5 4 
Segment 



TOP, bottom layer 



Measured 
Simulated 



18 17 16 



1 



Segment 



(d) 



Fig. 12. Measured occupancies together with the best predictions for the irradiation 
with 2°^T1 from the top for (a) 0.58 MeV, middle layer (b) 2.61 MeV, middle layer 
(c) 0.58 MeV, bottom layer (d) 2.61 MeV, bottom layer. 

ties. The background has to be well known, if lines are used that also occur 
in the background. In general, a source should be chosen which has lines with 
a reasonably high probability of containment in a segment and high enough 
energy not to be disturbed too much by the surrounding material. The selec- 
tion will depend on the setup and the detector. In the case considered here, 
the ^"^Co lines were found to be best suited. 



5 Conclusions 



A new method is presented to determine the axes orientation of (/(-segmented 
germanium detectors. The method relies on the ability to predict occupancy 
patterns caused by the transverse anisotropy of the drift of the relevant charge 
carriers. The axes orientation of a test detector was determined with a preci- 
sion of about 6°. The precision was limited by the knowledge about the source 
locations. In general, a precision of 5° should be attainable. This is compat- 
ible to the precision achieved with the well known scan method and is the 
precision needed for the usages of pulse shapes as currently envisioned. The 
new method can be used in any setup, also with multiple detectors, using 
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Fig. 13. Measured occupancies together with the best predictions for the irradiation 
with from the side for (a) 0.58 MeV, middle layer (b) 2.61 MeV, middle layer 

(c) 0.58 MeV, bottom layer (d) 2.61 MeV, bottom layer. 

any calibration source as long as the setup can be simulated with sufficient 
accuracy and background, if present, is properly taken into account. There- 
fore, it becomes unnecessary to perform time consuming scans to characterize 
detectors before they are integrated in a larger device. The method was de- 
veloped and tested for true-coaxial detectors. However, it can be applied to 
any detector with an appropriate segmentation scheme. The modifications for 
a 0-segmented closed-end cylindrical detector will be minimal. 

The method has systematic uncertainties depending on the precision of the 
simulation of the setup and of the transverse anisotropy in the detector crystal. 
The dependence of the patterns on the transverse anisotropy provides the 
potential to determine charge carrier mobilities. 



6 Acknowledgements 



We would like to thank the members of the MaGe Monte Carlo group for 
their kind support. 



17 



References 



[1] G. F. Knoll, Radiation Detection and Measurement third ed., (Wiley, New York 
1999). 

[2] J. Eberth and J. Simpson, Prog. Part. Nucl. Phys. 60 (2008) 283. 
[3] I-Y. Lee and J. Simpson, Nucl. Phys. News 20 Issue 4 (2010) 23. 



[4] I. Abt et al, GERDA: Letter of Intent, hep-ex/0404039 



[5] F. Recchia et al, Nucl. Instr. Meth. A 604 (2009) 555. 

[6] I. Abt et al, Nucl. Instr. Meth. A 583 (2007) 332. 

[7] I. Abt et al, Eur. Phys. J. A, 36 (2008) 139. 

[8] B. Bruyneel et al, Nucl. Instr. Meth. A 569 (2006) 764. 

[9] I. Abt et al, Eur. Phys. J. C 68 (2010) 609. 

[10] I. Abt et al, Nucl. Instr. Meth. A 577 (2007) 574. 

[11] S. R. Elliott et al, Nucl. Instr. Meth. A 558 (2006) 504. 

[12] I. Abt et al, Eur. Phys. J. C 52 (2007) 19. 

[13] I. Abt et al, Nucl. Instr. Meth. A 570 (2007) 479. 

[14] I. Abt et al, J. Instr 4 (2009) 11008. 

[15] I. Abt et al, Eur. Phys. J. Appl. Phys, 56 (2011) 10104. 

[16] M. Bosweh et al, IEEE Trans. Nucl. Sci. 58 (2011) 1212. 



18 



